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Abstract 


There has been an extensive research going on in the field of graphical simulation 
of manufacturing processes. Most often such simulations rely almost entirely on 
numerical methods. There is another, a symbolic, approach in which the values to 
variable entities are bound in the final steps, and the intermediate calculations are 
done in a symbolic fashion. 

In this work such a simulation of machine tooling process is done. For this pur- 
pose a symbolic manipulator has been designed and developed which manipulates 
the constituent equations in a symbolic manner, and then the output, which is the 
mathematical counterpart of the actual surface, is rendered on a graphics worksta- 
tion. Several examples have been taken and directly tested and the results were 
satisfactory. 

The user is provided with a textual menu which prompts him to enter specifications 
for a cutter, the kinematics of motion of the blank and the cutter and some rendering 
options. The program then produces the output which he can see on the terminal. 
This package allows a designer to view a machining process before it is actually done. 
This enables him to avoid pitfalls, if any, and take corrective actions. This also 
saves a lot of time, effort and money which are invariably associated with the actual 
machining. 

The present work finds applications in tool and cutter grinding, gear cutting, 
turbine blade cutting and many other similar machine tooling processes. 
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Chapter 1 
Introduction 


1.1 Geometric Modeling 

The modeling of the geometric and topological features comes under the field of geo- 
metric modeling. Geometric information is essential for describing products, because 
every machine functionality is realized one way or the other by its geometric features. 

There are basically three levels of geometric modeling based on the dimensionality 
of the representational domain. These levels are : 

• Wireframe modeling. 

• Surface modeling. 

• Solid modeling. 

Out of these three, solid modeling has the most explicit representation of the 
objects in the real world and implicitly embodies the lower levels of representation as 
well. To make the concept clearer, we give a simple example. 

We know that the equation 

explicitly represents a solid sphere centered at origin and with radius r. It is easy 
to see that the same equation implicitly contains the equation of the surface of the 
same sphere.( Just change the < to =.) This was a very primitive example. In general 
it is not necessarily true that any solid modeling approach can uniformly manipulate 
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cut it it's of low('r diincnsioiiality as well as solids. Formerly there was a notion that 
for prcxlnct shape description solid modeling framework is the best. However, it has 
b«‘n realiz<‘d that solid modeling is not the panacea for handling product shapes. 
For example, wireframe models are useful for analysis of mechanisms wherezis surface 
modeling is appropriate for designing of maehining applications. In fact, the present 
work relics on surface modeling extensively. 

1.1.1 Conjugate Surfaces 

In many situations, we encounter surfaces which are not independent but related to 
one another by some constraint. This is more so in the case of machining processes, 
where one surface is constantly moulding the other. In the present situation, the 
surfaces and Jlia related by what is known as the higher pair contact constraint 
also known as conjugate surface constraint. These are the surfaces which have a 
common normal at the point of contact. 

In cutter path simulation the two surfaces ]Ci and are conjugate to one another 
b<‘cause the machined surface is nothing but different instances of the cutter profile. 
'I'hi.s will h<‘coine clearer in subsequent chapters. In practice, however, there will be 
slight deviations because of the non-ideal nature of the cutter and blank materials. 

1.2 Literature Review 

One of the major areas of research in manufacturing is the simulation of machining 
processes. The two basic determinants of a machining process are the cutter geometry 
and the geometry of the relative motion between the cutter and the work piece. 
There are many kinematic chain models of relative motion in literature. Srinivasan 
and Ferreira (1990) have proposed a two level model for machining processes. The 
lower level models the kinematic chain and the higher level decomposes the relative 
motion into cutting and feed motion components and with each component a subset of 
kinematic chain from lower level is associated. They have shown that for a significant 
subset of all machining processes their approach is useful in feature definition and 
validation, process selection and kinematic planning at the machining center. 
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There ha.s been a lot of development in the area of geometric modeling, particulaxly 
solid modeling which has provided design and manufacturing engineers with elegant 
and precise definition of components or products that axe to be designed (Mortenson 
1985; Mantyla 1988; Voelcker and Requicha, 1977; Faux and Pratt, 1988). Based on 
these definitions, several research workers have developed methodologies to generate 
process plans for manufacturing the designed components (Chang and Wysle, 1985; 
Kanumury and Chang, 1991; Joshi and Chang, 1988; Perng et. al,1990; Shah, 1991; 
Srinivasan and Ferreira, 1990). Most of the process plans are based on the geometric 
features of a component and the capability of the manufacturing process to produce 
those features. A lot of research has been done regarding the physics and mechanics 
of several processes (Acherkan et. ah, 1973; Ghosh and Mallik, 1981; Ber et.al., 1988; 
Sen and Bhattacharya, 1969; Pande and Shah, 1980) but most of it is directed towards 
the area of gear manufacturing. 

Wang and Wang (1986 a, 1986 b) have defined swept volume using envelope 
theory approach. However they don’t model the kinematic structure of the machine 
tool explicitly. Yap (1988) has given some algorithms for computing 3-D tool paths 
for planar and non-planar machined surfaces. Dhande et al. (1990, 1992) have shown 
how the gtx)metry of a machined surface can be derived from the wire-EDM, which 
is a non-traditional manufacturing process. Drysdale and Jerard (1989) have done a 
discrete simulation of machining for sculptured surfaces. Blackmore and Len (1990) 
have proposed a swept differential equation approach to define the swept volume. 
Generalized sweep operations are discussed by K.Weiler and D.Melachlan (1990). 

The problem of conjugate geometry between a pair of moving as well as contacting 
surfaces is termed as higher pair contact problem in kinematics (Boltyanskii, 1964; 
Dudley and Darle, 1962; Chakraborthy and Dhande, 1977). These problems have 
been analyzed for design of cams and gears. Dhande and Chakraborthy, 1977, have 
given methods to define kinematics of relative motion as well as the geometry of 
contacting surfaces for 3-D cams. Some work has also been done concerning the 
conjugate geometry of planar and spatial gear mechanisms (Dudley, 1962; Dhande 
and Chakraborthy, 1977). 
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One of the major coniputational difficulties in formulating and analyzing the prob- 
lems of conjugate geometry is that for every given situation of a pair of contacting 
surfaces and their relative motion, one has to derive the biparametric equations of 
the surfaces every time. This can be improved by using the symbolic manipulation 
approach for the problems of conjugate geometry , especially for modeling different 
manufacturing processes. Reinholtz and Dhande et aL(1990) have described an ap- 
proach of how symbolic manipulation can be used for defining the machined surfaces. 
I hey have considered only a three degrees of freedom model to show how helically 
swept surfaces can be manufactured using alternative strategies. 

1.3 Objectives, Scope and Organization 

The objective of the present work is to graphically render the path taken by the cutter, 
given the parameters describing its surface and kinematics regarding the relative 
motion of the cutter and the blank with all the intermediate computations done in 
symbolic form. 

This report consists of five chapters. Chapter 1 briefly introduces the field of ge- 
onretric mcxkding and conjugate geometry followed by an extensive literature review. 

Chapter 2 explains the problem, drawing analogy from 2-D envelope theory. Then 
the symbolic model of the generic cutter and the symbolic algorithm to derive the 
cutter path (in other words, the machined surface) are presented. And finally the 
computational part of the algorithm is explained followed by some observations on 
the limitations of this (envelope theory) approach. 

Chapter 3 starts with the role of the symbolic computation in the present work 
and goes into the design and implementation aspects of the symbolic manipulator 
specifically designed for this work. At the end some examples are given to illustrate 
the capabilities of the symbolic manipulator. 

Chapter 4 is about the rendition of surfaces on a graphics workstation. It also 
discusses the interface provided to the user and some tips on how to use it. It 
concludes with with some examples which have been tested and run on a graphics 
workstation. 
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C'hai>t(T 5 conclude's the work with 


some technical details and suggestions for 


further extension of this work. 


Chapter 2 

Mathematical formulation 


2.1 Envelope Theory - 2D 

In this section we w'ill give a brief summary of the important results of the envelope 
theory. 

A family of curves can be specified by the implicit form equation 

f{x, y, a) = 0 (2-1) 

where x and y are the Cartesian coordinates of a generic point and a is a parameter 
(S<K* Fig.2.1). By giving a specific numerical value of a, say one can obtain the 
equation of a member of the family as 

/(x, y, a.) = 0 (2.2) 

It is aissumed that every member of the family is a C^-continuous (A -continuous 
curve is one which has first differential at all points on it) curve. Furthermore, it is 
assumed that successive members of the family intersect at a finite point. Each of 
these points of intersections, it is assumed, has a limiting position when the successive 
member curves are considered infinitesimally close and approaching a limit (Boltyan- 
skii, 1964). The point of intersection of a pair of infinitesimally close member curves 
can be obtained by solving the equations 


/(x, y, a.) = 0 


(2.3) 
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f{x, y, Qri + Aa) = 0 


(2.4) 


Afi being an infinitesimal increase in a. The limiting position of this point of 
inter.section ran be obtained by eliminating a from the following equations: 


fix, y, a) = 0 (2.5) 

Is = » 

All such linuting positions, in turn, define the envelope to the family of curves. 
In parametric form, this equation of a family of curves is given by 


/(u, a) = 0 


and a generic point P of this can be expressed as 


(2.7) 


p{u, q) - [x(u, a) y(u, q)] (2.8) 

wh<‘re Umin < « < timux Qmin < » < ocmax- In parametric form, the limiting 
points of intersections can be obtained by considering the collinearity of the vectors 

au da 


In other words, the equation of the envelope in parametric form can be obtained by 
eliminating the parameter ti between the following two equations (Chahraborty and 
Dhande, 1977): 


f(u, a) = 0 

(2.9) 

n-V = 0 

(2.10) 


where 


- dp ,dps 

V = and n • (^) = -1 - 
da ou 
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It will 1)0 in place to give a physical interpretation to Equation (2.9). If the 
paraniotx'r o represents time then V and n will be the velocity and normal vectors at 
any point. When a curve is being swept, every point on the curve is moving along 
the direction of its velocity. If this velocity vector is orthogonal to the normal vector 
of the curve at any point (the C^-continuity ensures that there would be atleast one 
such point) then that point of the curve is the point defining the envelope curve at 
that instant. All other points will lie within the swept area of the curve. 

I'he above description was for the envelope of a curve in a fixed frame of reference. 
Many tinit's it is required to find the envelope of a curve or a surface in another frame 
of reference which is moving with respect to the frame of reference in which the curve 
is moving. In chapter 3 we talk about one such case. 

2.1.1 Symbolic Model 

In any machining process there are three elements which describe the machining 
process completely (assuming ideal behavior of cutter and blank material). They are 


1. The surface of the cutting tool (denoted by throughout this work). 

2. The surface to be produced, denoted by ^2 *■^<1 

3. I'he relative motion between the above mentioned surfaces. 

To be able to produce the the surface E 2 properly, it is necessary to assume that 
the two surfaces Ei and E 2 are conjugate to each other. Or in other words, these 
two surfaces maintain a higher pair contact throughout the cutting motion. In the 
following sections we will give the symbolic models for the generic cutter (Ei) as well 
as the developed surface (E 2 )- But let us first recapitulate the symbols used in the 

derivations. 

2.1.2 Basic Concepts 

A generic model of a machine tool is essentially a representation of the kinematic 
structure of the machine tool ( Acherkan, 1973). The kinematic structure of a machine 
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tool describes the kinematics needed for giving the cutter its cutting motion, and the 
feed motions to the cutter and/or the blank. 

In the present model (Karunakara Poopathi, K.P, 1993), the cutter and the blank 
are considered to be free bodies in space and therefore their positions and orientations 
(w.r.t the coordinate axes of a fixed frame of reference) can be described by twelve 
functions (of time), six for each surface. The fixed frame of reference is denoted by 
So (Oo — Xo, Vo, Zo). A coordinate system Si (0i — Xi, Pi, Zi), is assumed 
to be attached to the cutter, and another coordinate system S 2 (O 2 — X 2 , Y 2 , Z 2 ), 
to the blank. The translation of the origins 0i and O 2 of the moving coordinate 
frames Si and S 2 respectively are (ai,6i,Ci) and (a 2 ,i» 2 ,C 2 ) whereas (^ 1 ,^ 1 , V’l) and 
denote the rotation of the reference frames around the three axes of the 
fixed coordinate frame. 

2.2 The Generic Cutter 

2.2.1 Introduction 

The geometry of a generic cutter can be modeled as an axisymmetric, biparametric 
surface. The axis of symmetry depends on whether the cutter is end mill type cutter, 
a side mill type cutter, a face mill type cutter or any combination of these basic types 
of cutter. 

In numerical control literature, a generic cutter is defined using seven parameters 
and the generatrix consists of two straight line segments forming the tip flank as well 
as the side flank and a circular fillet between these straight line segments (Chang, 
1989). In this work, an extension of this is used (Karunakara Poopathi, K.P, 1993). 
The generatrix curve is assumed to be a piecewise continuous curve consisting of three 
line segments and a circular fillet between the tip flank and the side flank lines. This 
generatrix curve can be defined by means of eight parameters, \iz.,d, r, e, f, a, b, h 
and hi. The cutter surface can be obtained by rotating the generatrix curve either 
along Zi axis or along the axis. Rotation about Zi axis produces various types 
of end mill cutters wherecis rotation around Xi axis generates a disc type cutting 
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tool geometry which will represent a side-and-face milling cutter or a grinding wheel. 
If the directrix is helical, then the geometry of a gear hob, or a tapping tool can 
be obtained. Thus by choosing appropriate values of the parameters describing the 
planar curve and the sweep motion, one can obtain a variety of conventional and 
non-traditional cutting tools. 

2.2.2 Description of the generic cutter 

The planar curve of the generic cutter surface consists of three straight line segments 
and a circular arc segment. For the case of end mill type cutter, the first and second 
line segments will generate the tip and side flanks of the cutter while the third straight 
line segment will generally represent the shank of the cutter. The circular segment 
will generate the corner radius of the end mill. For this purpose, a coordinate system 
5c (Oc - n, ^c) is considered to be attached to this cutter profile. This planar 
profile of the cutter can be represented by means of a parameter u which is the 
profile length of any point on the profile with respect to origin Oc- The cutter surface 
can be obtained by sweeping the this generatrix profile along a directrix curve 
which can be circular, helical etc. This sweep motion is represented by by means 
of a homogeneous transformation matrix [ \M\ whose elements will be functions of 
parameter u. 


OcP : Zc = Xc tan a 

(2.11) 

PQ:{xc-ef-\-{zc-f) = r^ 

(2.12) 

QS : Xc = Zc tan tan a tan 6 + - 

(2.13) 

d d 

ST : Xc = h tan b — — tan a tan 6 + — 

(2.14) 


where 


OcP - First straight line segment, 
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PQ - The circular segment with center at A, 

QS - Second straight line segment, 

ST - Third straight line segment, 

R - Point of intersection of the segment OcP and the third segment QS, 
d - Diameter of the cutter measured at R, 

r - Radius of the circular segment. This is positive for a convex arc and 
negative for a concave arc, 

a - Angle of segment OcP measured from Xc axis. 

— X X 

b - Angle of segment QS measured from Zc axis. 


— X 


2 


<b< 


IT 

2 ’ 


e - X coordinate of point A in the coordinate system Sc, 
f - Z coordinate of point A in the coordinate system Sc, 
h - Nonparallel height of the cutter, 
hi - Parallel height of the cutter, 

is.Jilletj'adius - Flag which is to be set to TRUE when the circular 
segment PQ is tangent to both the straight line segments OcP and QS and FALSE 
otherwise, 

u - Parameter which is the profile length of any generic point P on the 
profile measured from Oc, 

V - Parameter which describes the sweep of the two dimensional profile, 

[ \M\ - Homogeneous transformation matrix which describes the sweep of 
the two dimensional profile, 

51 - Length of the first straight line segment OcP, 

52 - Length of the circular arc segment PQ, 

53 - Length of the second straight line segment QS, 

54 - Length of the third straight line segment ST, 
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- Position vector of a generic point P on the two dimensional profile 
of the cutter in Sc coordinate system in terms of u and 

- Position vector of a generic point P on the surface of the cutter in 
Si coordinate system in terms of u and v. 

The parametric equation of the generatrix curve with respect to the 

coordinate system Sc{Oc — Xc,Yc.,Zc) can now be defined as : 


V(u,u) = ‘p(u)-[lM] (2.15) 

where [ JM] is the matrix describing the sweep motion. 

2.2.3 The Algorithm 

In this section we will give the algorithm used to derive the biparametric equation 
of the surface from the nine parameters of the cutter profile described earlier 
(Karunakara Poopathi,K.P, 1993). It has been implemented in the symbolic manip- 
ulator designed in this work. It may be noted that all the parameters are not needed 
in all cases. The algorithm may be stated as follows ; 

1. Calculate e and / if is^filletjradius = TRUE. 


d ( cos a — sin h 
2 cos(a-h' 


h) 


(2.16) 


/ 


r + e sin o 
cos a 


(2.17) 


2 . 


Calculate point P. 


= (e + / tan a) cos' a- 


■\/[{e + f tan a) cos^ a]* — (e^ + /^ — r* cos^ a (2.18) 


= xf tan a 


(2.19) 


3. 


Calculate point Q. 



toi 
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= {/ + [e -f ^(tanatan 6 — 1)] tan &} cos^ 6 + 

(tan a tan — 1 )] tan b) cos^ 6] — {/^ + [e + ^(tan a tan b — 1)]^ — r*} cos^ b 

£4 

( 2 . 20 ) 



= zf tan b — -(tan a tan 6—1) 

£ 

(2.21) 

4 . 

Calculate h if is.filletjradius = TRUE and 6=0. 



h = zf 

(2.22) 

5 . 

Calculate point S. 



s 

= h tan ^ — 2 ^ ^ 2 

( 2 . 23 ) 


II 

( 2 . 24 ) 

6. 

Calculate 5 i, 52, 53 and 54. 



5 l = 

( 2 . 25 ) 


52 = r / arctan ( -7^ | arctan [ ) 1 

i \x 7 -ej \^c-^J) 

( 2 . 26 ) 


= v/(^f - 

( 2 . 27 ) 


S4 = hi 

( 2 . 28 ) 


7 . Define the parametric equation 


y{u) = [Xc(w) yc{u) Zc{u) 1 ] 
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where, 

for the interval 0 < u < Si 
Xc{u) = u cos a 

Vein) = 0 ( 2 . 29 ) 

2c(u) = u sin a 

for the interval Si < u < 

Xc{u) = e + abs{r) cos larctan 

ydu) = 0 ( 2 . 30 ) 

/ \ r T / \ • f /sisina — /\ /u — si\l 

zJu) = f + abs(r) sin < arctan I | + | I > 

( \ Si cos a — ej \ r J ) 

for the interval Si + sj < u < 5 i + S2 + 53, 

Xc(«) = [u - (si -f S2 + 53)] sin b + x^ 

VcH = 0 ( 2 . 31 ) 

Zc(«) = [tl - (Si -f 32 + 33 )] COS 6 + zf 

and for the interval 31+S2+33 < « < Si+S2+'S3+34, 
Xc{u) = xf 

ydu) = 0 ( 2 . 32 ) 

zdu) = [u - (si + 32 + 33 )] + zf . 

8 . Define the sweep matrix [ \M] 

9 . Define ^^t£,u) 


Si sina — / 
Si cos a — e 


IX — Si 


y(u,u)= ''p{u)-[l^] 


( 2 . 33 ) 
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We give some examples of descriptions of the generic cutter. The input parame- 
ters, the symbolic output of the symbolic manipulator and the corresponding profile 
of the cutter are shown. The data are manipulated so that they are in more readable 
form. 


2.3 The Machined Surface 

In the earlier sections we saw how a generic cutter is described and how the various 
motion parameters are represented in the symbolic model. In this section we will 
give the necessary equations to synthesize the symbolic equations for the machined 
surface Yli- 

As the equations of the generic cutter and the kinematics are known, it is possible^ 
to derive the symbolic equation(s) for the path of the cutter, which is nothing but 
the required machined surface X)2- The cutter surface is expressed as a biparametric 
surface in u and v. As the two surfaces are assumed to be ideal and in higher pair 
contact, at the point of contact P we have, 

;mi = VI (2M) 

and 

'n[°M] = ± 'n[5M] (2.35) 

where, 

- Point P with reference to the cutter coordinate system Si- 

- Point P with reference to the blank coordinate system S2. 

^rip^i - Unit normal vector of the cutter surface point P with 

reference to coordinate system S^. 

- Unit normal vector of the blank surface ^^2 point P with 
reference to coordinate system S2- 


^We will see in subsequent discussion that it may not be easy in all cases. 
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Profile of the Cutter 


Is arc 

tangential? 

■ 

■ 

a 

b 

mn 

e 

f 

■ 

Yes 

25.0 



0.0 

57.5 

0.0 




Input Data to Symbolic Manipulator 


Range of u 

Xc 

Yc 


- 

- 

- 

- 

0.0< u < 19.635 

12.5cos[0.08u- 1.571] 

0.0 

12.5sin[0.08u- 1.571] 

- 

- 

- 


19.635 < ti < 77.135 

12.5 

0.0 

u - 7.135 


Output from Symbolic Manipulator 


Figure 2.3: Ball Nose End Mill Cutter 
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Profile of the Cutter 


Is arc 
tangential? 

d 

1 

a 

■ 

hi 

e 

f 

h 

Yes 

25.0 



0.0 

55.843 

8.5 

14.157 

5.561 


Input Data to Symbolic Manipulator 


Range of u 

Xc 


Zc 

0.0< u < 16.021 

0.707u 

■ilil 

0.707u 

16.021< « < 19.162 

4 cos[0.25(«- 16.021) 
-0. 785] -f 0.508 

nn 

4sin[0.25(u - 16.021) 
-0.7851+14.157 

- 

- 

- 

- 

19.162 < « < 75.006 

12.5 

0.0 

u — 5.006 


Output from Symbolic Manipulator 


Figure 2.4: 90 deg. Tool 
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Profile of the Cutter 


Is arc 

tangential? 

d 

■ 

a 

■ 

hi 

e 

f 

h 

Yes 

25.0 

4.0 

10.0 

5.0 

50.0 

8.782 

5.61 

20.0 


Input Data to Symbolic Manipulator 


Rajige of u 


Yc 

Zc 

0.0< u < 9.624 

0.985U 

0.0 

0.174U 

9.624< u < 14.860 

4 cos[0.25(u — 9.264) 
-1.396]+8.783 

0.0 

4 cos[0.25(u — 9.264) 
-1.396]+5.61 

14.86< u < 29.654 

0.087(u-29.654)+14.057 

0.0 

0.996(u-29.654)+5.61 

29.654 <u< 79.654 

14.057 

0.0 

ti — 9.654 


Output from Symbolic Manipulator 


Figure 2.5: Side and Bottom Angle Cutter 
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Profile of the Cutter 


Is arc 

tangential? 

II 

■ 

a 

■ 

■■ 

e 

f 

■1 

No 

25.0 

6.0 

0.0 

0.0 

61.528 

8.5 

4.0 

8.472 


Input Data to Symbolic Manipulator 


Range of u 



Zc 

0.0< u < 4.028 

u 

0.0 

0.0 

4.028< u < 23.545 

6cos[0.167(tt-4.028) 

-2.412]-h8.5 

0.0 

6cos[0.167(u — 4.028) 
-2.412]-t-4.0 

- 

- 

- 

- 

23.545 < u < 85.073 

12.5 

0.0 

u - 15.073 


Output from Symbolic Manipulator 


Figure 2.6: Non-Fillet Type End Mill 
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Profile of the Cutter 


Is arc 

tangential? 

d 

■ 

a 

b 

|^H| 

e 

f 

■ 

Yes 

5.181 

4.0 

0.0 



0.0 

4.0 

50.0 


Input Data to Symbolic Manipulator 


Range of u 




- 

- 

- 

- 

0.0< u < 4.608 

4cos[0.025ti- 1.571] 

0.0 

4cos[0.025u-1.571] 

4.608< « < 56.809 

-0.837(u - 56.809) + 25.0 

0.0 

0.548(u-56.809)+50.0 

56.809 < u < 76.809 

25.0 

0.0 

u - 6.809 


Output from Symbolic Manipulator 


Figure 2.7: Conical Cutter 
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[ \M] - Homogeneous transformation matrix for transforming a vector 
from the coordinate system Si to Sq. 

[ 2 ^/] - Homogeneous transformation matrix for transforming a vector 
from the coordinate system to So- 

The surface ^2 is the boundary of the swept volume of the cutter surface J2i, 
which is nothing but the envelope to the various instances of the cutter surface. We 
observe that at any given instance, only a single profile is cutting the blank surface^. 
We also observe that this particular profile (out of potentially infinite profiles) is 
the one where the common normal (the surface are assumed to be conjugate) of the 
surfaces and YI 2 is orthogonal to the relative velocity vector between the cutter 
and the blank (Chakraborthy and Dhande, 1977). This is what is technically known 
as condition of cutting or the condition of contact or the condition of envelope. 
Mathematically, it can be expressed as. 


* ^kp,i2 — 0 (2.36) 

and 

“Sp., ■ “VV,,2 = 0 (2.37) 

where, 

“«p,i - Unit normal vector of the cutter surface 
at point P with reference to coordinate system ^o- 

®np ,2 - Unit normal vector of the blank surface E 2 
at point P with reference to coordinate system 5o- 

°Vpi 2 - Relative velocity between the tool and blank with reference to coordinate 

system Sq. 

The condition of contact, being a dot-product between two vectors, is a scalar 
equation relating parameters u, t; and t. Thus, at any given instant of time t, the 
condition of contact reduces to a relation between u and v. For any value of u within 


^Later, a case will be seen, where this is not the case. But for most cases this holds true. 
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its fecisible range, one can find the corresponding value of the parameter v. Thus, 
one can locate points of surface JZi which define the surface '^2- The corresponding 
points of the surface XI 2 ^^.n be obtained using the equation, 

V= ^P,\M][lM] (2.38) 

2.3.1 The Algorithm 

In this section, we give the algorithm to find the machined surface J^2- h has been 
programmed in the symbolic manipulator developed in this work. It can be stated as 
follows : 

1 . Define the position vector of a generic point P on the cutter surface 

as a biparametric surface in u and v such that 


^p(u,v) = [xi(«,v) yi{u,v) zi(u,v) 1] (2.39) 

2. Define the twelve parameters of motion Cl, fej, ci, 61, 0i, 02 , I> 2 » C 2 , 

62, <l>2 o.nd in terms of t. 

3. Calculate the matrices [ jM] and [ ^M] as follows : 


= (2.40) 

where. 


[ w = 

[W = 



■ 1 0 0 0 ' 

0 cos 9i sin 0 

0 — sin cos 9i 0 

. 0 0 0 1 . 

’ cos (f>i 0 — sin <j>i 0 ' 
0 10 0 

sin 0 cos (^1 0 

.0 0 0 1 . 

cosV»i sin^’i 0 0 ' 

— sinV’i cosV’i 0 0 

0 0 10 

0 0 0 1 
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1 0 0 0 ' 

0 10 0 

0 0 10 

flj ^ ^1 1 _ 

Matrix [ °M] is the time-differential of matrix [ jM]. 
4. Calculate the matrices [ jAf] and [ as follows : 



where, 




[M 


[X] 


[X] 


10 0 0 

0 cos $2 sin $2 0 

0 — sin $2 cos $2 0 

0 0 0 1 

cos (j >2 0 — sin (l >2 0 

0 10 0 

sin <{>2 0 cos (j >2 0 

0 0 0 1 


cos ^2 

— sin ^>2 
0 
0 


[IT] 


10 0 0 
0 10 0 
0 0 10 
02 &2 ^2 1 


sin ^2 0 0 
cos ^2 0 0 
0 1 0 
0 0 1 


EM) = gA/]-- 

Matrix [ “Mj is the time-differential of matrix [ 5^]* 

5. Derive the position vector ^p. 

= [X2 J/2 22 1] 

v= 'PSm&M] 

6. Derive the normal vector ^Np,i. 


(2.41) 


(2.42) 
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where, 


'«» ■ £ « 


d^p 

dv 


( 2 . 43 ) 


9V . 

— ond ^ 

UU OV 


are the partial derivatives of with respect to u and v respectively. Therefore, 


d^p _ dyi dzx , 

du du du du 


d^p 

dv 


When the vectors 


dyi dzj 
dv dv dv 


d^p , d^p 
and 

du dv 


i 

j 

k 

dx% 

§20. 

dz\ 


9yi 


dv 

dv 

dv 


are expressed in cartesian coordinates, the normal vector ^Np,i is defined as 

= 

The same can be redefined when the vectors 

d^p j d^p 
du av 

are expressed in homogeneous coordinates as : 

* j 


^Np,i = 


dxi 

§1LL 

dzi 


& 


dv 

dv 

dv 

0 

0 

0 

1 vector® 



np,i 


'Np,, 


( 2 . 44 ) 


^Actually this step may be skipped. 
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where 

is the magnitude of 

8. Calculate the unit normal vector'* 

np,i 

in the global coordinate system. 

= 'np.x [ W (2.45) 

9. Calculate the relative velocity vector °Vp,i 2 - 

= “Vp,2 - ^Vp,r (2.46) 

where, 

°Vp.i = ^p[°M]and ^Vp,2 = ^p[lM] 

10. Reduce symbolically the condition of contact 

■ *^Vj=>,i2 = 0 (2-47) 

11. obtained in step 5 will be a function of u,v and t. Use the condition of 
contact obtained in step 10 to eliminate either u or v. This will yield the surface Yli- 

■*This will not be unit normal if the earlier step is skipped. But that makes no difference as we 
will observe later. 
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2.4 Computational Model 

As observed earlier, it may not always be possible to symbolically eliminate a variable 
between equations (2.42) and (2.47). Infact, even for the most simple kinematics 
the equations become very complex. In such circumstances we revert to numerical 
methods. 

As explained earlier, the generic cutter is modeled as a biparametric , axisymmet- 
ric surface. Equation (2.42) gives us the cutter surface at any time instant t in the 
coordinate reference frame 82 - Hence it is a function of u, v and i. Similarly, equa- 
tion (2.47) gives the condition of contact which the cutter surface must satisfy at any 
time. This condition of contact depends upon the cutter geometry and the relative 
kinematics between the cutter and blank surfaces. A point worth noticing is that the 
condition of contact is independent of the coordinate reference frame chosen. This 
does not mean that the equation for condition of contact is same for all coordinate 
systems, but only that the solution value of t; is same in all coordinate systems. This 
is because neither of the transformations, viz. translation or rotation changes the 
angles between the vectors and °Vp,i 2 • This condition of contact is also a 

function of u, v and t. Now , to get the surface '£,2 we have to eliminate either of u 
or V. It can be easily observed that it would be most convenient to eliminate v. We 
also make the following observations : 

• For a given time instant t' and parametric value u' of u, there are two values of 
V, say Vi and V 2 which satisfy the condition of contact and such that 

Ul = V2 + TT 

• For such ui, V 2 and t' and for all u the points ^p{u,vi,t') and ‘^p(u,V 2 ,f) are 
on the surface X) 2 - 

All the corresponding points of the surface X )2 then calculated and stored in 
an array and rendered later. 
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2.5 Some Observations 

Till now we were tacitly assuming that there will be only a single profile of the cutter 
which will be orthogonal to the relative velocity vector. This is not the Ccise. Even 
if this were the case, there are other situations in which this algorithm fails. The 
failure does not imply any fault in the symbolic equations derived as such, but that 
such equations cannot always be rendered automatically. In such places the user 
intervention becomes necessary. And the user needs to be a C - programmer well 
versed in the use of graphics package employed (in this case - STARBASE). In the 
following discussion we give some cases where user intervention becomes necessary : 

1 ) Not all types of kinematics are allowed : There are many simple cases of 
motion of the cutter and blank where the present approach fails. One such case is the 
drilling motion. Assuming that the blank is stationary, consider the case when the 
cutter is moving along negative Z direction. Such a motion, in practice, will produce 
a hole in the blank. But in the present case we observe that there are not two but 
many (potentially infinite) profiles satisfying the condition of contact and thereby 
becoming eligible for inclusion in the machined surface. Even if we somehow show all 
the profiles the problem remains that how to connect each profile at a particular time 
instant with that of the next instant. And there is the problem of self intersection 
of the cutter because the previously cut portion of the blank by the tip of the cutter 
is cut later by the tail piece of the cutter. In such cases it becomes difficult (even in 
adhoc fashion) to determine which portion of the blank is to be retained and which 
removed. 

2) Crossing of the cutter path is not allowed : Even though ,logically, this 
is just another instance of the previous case, it deserves special attention. In cases of 
this type, everything goes on fine until ,suddenly, the cutter starts cutting the portion 
cut earlier. As earlier, here too it is not obvious what the final shape of the machined 
surface should be. Though there are some algorithms which calculate surface-surface 
intersection, it is not clear how much they will be useful in the present case. 

3) The cutter has to be axisymetric ; Throughout the work, it has been 
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assumed that the cutter is axisymmetric. Infact, this fact hais been used in the 
symbolic definition of the cutter, where only one half of the profile is defined and 
the rest of the cutter is developed by rotating this profile around the Z axis. But in 
practice this may not be the case, though even for most non-axisymmetric cutters we 
can argue that they are generally rotating, hence effectively forming an axisymmetric 
cutter surface. It should be noted that a cutter, such as a helically swept one, can 
be defined and used, but the rendition of final machined surface , is not easy. The 
problem here is that a helical cutter may screw its way into the blank and such cases 
are very tough to handle. 

4) The cutter motion is assumed continuous : Previously we said that the 
surface is obtained by joining different instances of cutter profile at different time 
instants such that the profile is orthogonal to the relative velocity vector. But in 
practice what generally happens is that a portion of the blank is machined and then 
the cutter is shifted to another portion and the machining started again. In such cases 
the cutter path (or the machined surface) is a set of unconnected surfaces instead of 
a single, continuous surface. How are such cases handled? 

In this work such cases are ignored. Such cases can be handled using 
interval functions. But we see that such cases can also be rendered by taiing different 
runs for each interval of the function. Incidentally, this is the flaw (the symbolic 
manipulator developed in this work does allow definitions of interval functions. But 
how to interpret them is left to the programmer®.) of the symbolic manipulator rather 
than that of the approach. 

5) The blank is assumed to be present everywhere : One of the assumptions 
which was made is that the cutter never runs out of blank, i.e the cutter is moving in 
a blank which practically extends to infinity in all directions. Obviously this is not 
the case in general. In many cases it becomes important to have an initial shape for 
the blank (All blanks have some initial shape) because in such cases we are interested 
in the piece that is left rather than the groove which is cut by the cutter. One such 
example is that of a gear. Though it is possible to visualize what the gear would look 

®For more information, kindly refer to the appendix-A given at the end of the report 
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like, it is taxing on the user’s imagination. 

And also in many situations we would like to have an initial shape for the 
blank just to see how many different workpieces could be made out of the given piece 
of material. If such a mechanism is provided, then the user could do the machining 
with different initial positions and orientations. One can then augment the basic 
simulator with one of the numerous algorithms for optimally fitting the workpieces 
into the given blank piece thereby minimizing the wastage of blank material. This is 
one of the most important criteria for the tool designer. 



Chapter 3 

Symbolic Manipulator 


3.1 Role of symbolic manipulator 

In this section we will discuss what a symbolic manipulator is, and what it does. After 
the initial history, we proceed to explain why at all it is needed and its advantages 
over conventional computing. 

3.1.1 Symbolic Manipulation — What is it? 

Solving problems with symbols substituted for actual entities is known as symbolic 
manipulation. In the present context, it means doing computation on variables which 
may or may not be bound to some particular numerical values. 

3.1.2 History and development 

Human species is unique in this respect that it is the only one to posess this capability 
to a substantial degree. Computers, which are replacing man in almost every field, 
have traditionally been used for numerical computation so much so that they have 
come to be known as number crunchers. But recently, with the advances in hard- 
ware and software technologies, it has been possible to use computers in symbolic 
computation. The idea, though, is not new. It dates back to 1840’s, when Lady Ada 
Lovelace^ first suggested it. 

^Lady Ada Lovelace was a patron of Charles Babbage, who is credited with the development of 
world’s first computer 
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Initially the symbolic manipulator packages were written for specific fields. But 
with the breakthrough in hardware technologies, we have a choice among a wide 
spectrum of computers, from P.Cs to mainframes, with computational speeds touching 
hundreds of Mflops(Mflops stands for Million floating point operations per second) 
and memories, correspondingly, running into several Giga bytes (1 Giga byte = 2^ 
bytes). With this, it has now become possible to develop packages which are general in 
nature and work on various platforms. MACSYMA and MathCAD axe two excellent 
examples. 

The symbolic manipulation systems are capable of performing a wide variety of 
operations. They can do not only polynomial, but also trigonometric and matrix 
algebra. They can solve a system of equations and find roots of an expression. They 
can also solve problems involving differentiation and integration. Some packages like 
MACSYMA can also find numerical solutions whenever needed. 

3.1.3 Why a symbolic manipulator? 

Today we are living in an age where graphical interfaces and visual computing are 
becoming increasingly common. The main reason behind their success is that they 
give more insight into the problem than any other mode of depicting information. 
One of the goals of any computation is conveying maximum information in minimum 
period of time. An algebraic equation, for example, conveys more information to a 
scientist or mathematician than the enumeration of many points of that equation. 

But this is not all there is to symbolic manipulator. Symbolic manipulator main- 
tains the generality of computation. Therefore, it is particularly well suited for devel- 
oping generic algorithms. Secondly, the precision of numerical systems is hardware 
dependent. The rounding and truncation may produce a cumulative effect and pro- 
duce erroneous results. With symbolic computation, the rounding and truncation 
can be deferred till the final steps, thereby minimizing the errors. Symbolic ma- 
nipulators can deal with zero, infinity and other such exceptional conditions more 
elegantly. Many times it is possible to use the intermediate computations for entirely 
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different purposes^. If it is desired to store the output, the storage overhead for nu- 
merical result is many times more than the corresponding symbolic one. Moreover, 
the numerical results are inextricably interwoven with the idiosyncracies of the data 
structures and numerical methods used, which is not the case with symbolic results. 

3.2 Software Design and Implementation 

The symbolic manipulator is the heart of the present work. It would, therefore be 
useful to study its inner workings. In this section we will be doing just that. 

Conceptually the symbolic manipulator consists of the following stages : 

• Parsing 

• Symbol table construction 

• Symbol substitution and function evaluation 

• Expression minimization 

It should be noted that the first two stages axe done in sequence for all the dec- 
larations and la^t two of the above steps are for each individual declaration as would 
become clearer later. 

Fig (3.1) gives a pictorial representation of the basic relationships among different 
modules. The actual relationships are however very complex. 

Parsing: In parsing stage the input is read one token at a time. For this purpose 
the parser calls the token recognizing routine repetitively and tries to fit the input into 
a predefined format called the grammar of the symbolic manipulator. The grammar 
defines what kind of input is expected and accepted by the symbolic manipulator.^ 
Most of the grammar resembles closely the C-grammar, so that there will not be much 
problem in using it. 

The main constructs recognized by the symbolic manipulator are: 

^One such case is discussed in chapter 4 

formal specification of the symbolic manipulator grammar is given in Appendix A 
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Figure 3.1: The Symbolic Manipulator 
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• Identifier declaration and initialization 

• Matrix declaration and initialization 

• Function declaration 

• Function call 

The initialization for functions and identifiers is done by a single expression and 
for matrices it is a comma separated list of expressions. An expression may contain, 
besides identifiers and constants, matrix and function names as well as function calls 
all connected by allowed operators. 

The parser constructs a parse tree in memory by using a dynamic array of tokens. 
This array is called Table. Each token contains information about the present entity 
stored there and the indices to left and right subtrees. The structure Token is defined 
as follows: 

struct Token{ 
double value; 
char string [20] ; 
int type; 
int ileft; 
int iright; 
int valid; 
int op; 

}; 


The field string holds the name of the token. The value field holds the numerical 
value if it is a constant. Field type holds the the type of the token, i.e. whether it is 
an identifier ,a constant or a matrix and so on. The fields ileft and iright are indices 
into the array where the left and the right subtrees may be found. Field valid was 
provided for the purpose of garbage collection. The last field op is the ASCII value 
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of the operators like +,-?/ so on and predefined values for operators like sin, cos 
etc. 

The root of the parse tree is left in the global variable out. Once the parsing is 
over, the tree is copied from Table to a similar array safe. The Table is then freed. 
The reason for this movement of the tree is that during computation, each individual 
declaration is brought into Table and manipulated. That leaves a lot of Tokens which 
are no longer used. To collect this garbage, the most effective way is to free the whole 
Table, which is possible only if the rest of the declarations axe stored in a safe place. 

Symbol Table Construction: Once the parse tree is in safe, a routine enter is 
called which, for each symbol encountered, makes an entry in the symbol table. In 
this process the latest definition overrides the previous one. It is also worth mentioning 
here that there should be no conflict of names even among different types of identifiers, 
e.g. between a matrix name and a function name. 

The symbol table is an array of structures called SYMJTAB which looks like: 


struct SYM_TAB{ 
chzo: ident [20] ; 
int defined; 
int type; 
int dep_list; 
int order; 
int definition; 
} 


The field ident stores the name of the token, if any. Field defined is non-zero 
if the corresponding symbol is defined and zero otherwise. The type of the token is 
given by type. For functions, depJist gives the variables on which the function might 
depend. The field order is exclusively for matrices and gives its order(through the 
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index into the safe) and definition is the index into the array where the definition 
of the current symbol may be found. 

The symbol table not only provides a convenient way to access information, but 
also associates a unique key (the slot number) with each symbol. This key can be 
used to develop more efficient algorithms, explained later. 

Substitution and Evaluation: After the symbol table is constructed the tree is 
again traversed and each declaration is copied into the Table and processed, one at a 
time. Once a declaration is in Table, its symbols are substituted by their correspond- 
ing definitions accessed from the safe through the symbol table. The substitution 
is done recursively, i.e. a symbol is first substituted fully before proceeding fur- 
ther. This unfortunately has a serious drawback (though it in no way affects the 
performance for the present purpose) that it reduces the capabilities of the symbolic 
manipulator somewhat. Later we discuss how this can be remedied with minimal 
changes in code. The problem is dealt with in chapterb. Here it would suffice to say 
that by cutting on this capability we have customized the symbolic manipulator to 
the present need. 

In the evaluation phase a function evaluate is invoked with current declaration 
(that which is in the Table at present) as its argument. This function as its name 
suggests, evaluates the function calls embedded in the current declaration. 

E.g. suppose that the current declaration is: 


Z * differentiate(x){i''2.0;}* differentiate(y){y''3.0;}; 


then evaluate would further invoke the differentiate routine with arguments x 
and ° (actually the indices into the Table ) followed by y and t/® ° and after that 
return. 

Minimization: After evaluation, the next step is minimization. There are several 
kinds of minimization performed. The major ones are: 

• Cancellation of terms with opposite(same) terms in addition(subtraction). 



42 


• Caiirellation of identical terms in division. 

• Cancellation of nested inverse operators like sin(arcsin(x + • - •)) i® reduced to 

X + . . .. 

• Cancellation of reverse operators in multiplication, like, sin(x) ♦ csc(x) is re- 
duced to 1.0. 

and so on. 

For this purpose, the operators are divided into the following groups: 

• The addition/subtraction group. 

• The multiplication/division group. 

• The trigonometric group. 

The trigonometric group triggers special treatment functions for trigonometric 
expressions and after that they also come into one of the first two groups. So, in 
effect, there are only two major groups. 

In the addition/subtraction group all expressions are treated as consisting of to- 
kens connected by operators binary -I- or — or unary or Sub-expressions 
consisting of other operators like *, / etc. are also 2 issumed to be a single albeit 
compound token. 

For this purpose a recursive function enlist is developed which is passed the root 
of a subtree. Depending upon the current operator at the root of the subtree, enlist 
decides which of the three groups the sub-expression falls into. Thus deciding on the 
group, it recursively calls itself on the left and right subtrees. 

At the bottom level of such call, i.e. when enlist reaches the leaves, enlist enters 
information about each of the leaves in an array of structures, called listnodes, A 
listnode looks like: 
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int key; 
int index; 
int valid; 
int direct; 
} 


The field key associates a key to each token which can be effectively used to sort 
them in case more efficient algorithms for minimization are used. The field index is 
the index into the Table where the token can be found. Field valid field is used in 
minimization to mark the cancelled tokens as invalid. The last field, direct, was used 
in the earlier stages of implementation and is no longer used. 

For each sub-expression we make two lists depending upon the group it belongs 
to. The addition/subtraction group has a plus- list and a minus-list. Similarly, the 
multiplication/division group has a num-list and den-list. These lists are stored in the 
same global array of structures listnodes. The array is an array of listnodes and is 
called list. The information regarding where each list starts and ends is passed across 
functions through pointer to structure listhead. A listhead consists of the following 
fields: 


struct listhead-C 
int 11; 
int 12; 
int type; 

} 


Field l\ is an index into array list and indicates where the plus-list or num-list 
starts. In a similar way 12 indicates where the other list starts. The end of the 
currently active (at any given time there are more than one set of lists except of 


44 


course when enlist is at the top level) minus or den list is stored in a global variable. 
Finally the field type is the indicator of the group the token belongs to. 

The bottom level enlist passes the pointer to this structure listhead to the routine 
minimize which does many minimizations and then enlist makes a new listhead, 
fills up the minimized information and returns it to the caller function. The caller 
function then extracts information from its left and right list heads merges them 
keeping in mind what types are the left and right listheads and proceeds further 
similar to the lower level enlist call. This goes on till the top level is reached and 
the minimization is over (note that the minimization is over ais far as this symbolic 
manipulator is concerned, it may not be, however, complete). Then the current 
minimized declaration is copied back back into the safe and the symbol table entry 
is changed to reflect the changes. 

Here it may be observed that a lot of time and memory may be saved using 
forward declaration of the symbols used, i.e. defining a symbol before it is used. 
However, this is not imperative. 

3.3 Some Pitfalls 

It is well known that any symbolic manipulator makes excessive demands on the 
system memory, as the intermediate expressions can become arbitrarily long. The 
symbolic manipulator developed in this work is no exception to this, though an at- 
tempt is made to reduce it by cutting on the generality of the symbolic manipulator. 

The design of the symbolic manipulator presented in the previous section is quite 
efficient with regards memory and run time^. But there are some cases, which had 
to be handled separately. 

Except for the functions inverse and taninv, all other functions are robust. Ro- 
bust in the sense that they can be used wherever it logically makes sense. Finding 
the cross product of two matrices makes sense, but crossing two functions doesn t, 

for example. 

‘‘For further improvements, refer to Chapters 
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The reason for inverse being fragile^ is that computing the inverse of a matrix 
symbolically is heavily taxing on the memory. The algorithm inverse, therefore, 
periodically flushes the Table in order to collect garbage. Hence, if the function call 
inverse occurs as on of the terms in an expression, then during the evaluation of 
inverse the rest of the expression is also flushed and lost. That is the reason why 
matrix inversions need to be done separately. 

The case of the function ianinv is different, taninv is similar to the standard C- 
library function, atan2, which takes two arguments, the Y and the X coordinate values 
and returns the corresponding angle in the interval [— 7r,+7r]. The need of taninv was 
realized in the final stages of the work, and hence only a minimal functionality was 
provided. The peculiarity is that taninv expects both its arguments to be computable 
at compile time. This however makes no difference present context, as they are indeed 
computable wherever it is needed in the present work. 

Another source of bugs might be the function declaration construct of the symbolic 
manipulator. In present application, nothing is to be gained by using it. It was 
provided for future extensions. As any usage of this construct can be modeled by 
other constructs, it was seldom tested thoroughly. 

3.4 Some Examples 

We illustrate the workings of the symbolic manipulator with some examples. These 
examples serve to show only some salient features of the symbolic manipulator and 
by no means are exhaustive. 

/♦INPUTl*/ 

y » differentiate(x){x''x + log(x);}- 2.0*x; 

z - (cos(acos(y)) + sin(y)*cosec(y))*tan(t)/(sec(u)*tan(t)*cos(u)) ; 


’^Robust and Fragile are terms borrowed from latex manual 
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matrix test [2 . 0 , 2 . 0] *{p , log (p) , 2 . 3. 4 . 9} ; 


testinv » inverse (){test ;}; 


/♦0UTPUT1*/ 

y- ((x)~(x))* (1. 000000+(1.000000)*(log(x) )) + (!. 000000)/Cx)- 
(2.000000)*(x) ; 

z«((x)‘(x))*(l.000000+(1.000000)*(log(i) ))+(!. 000000 
)/(x)+1.000000-(2.000000)*(x) ; 

test [2 . 000000 , 2 . 000000] {p , log (p) , 2 . 300000 ,4.900000 } ; 

test inv=mult [2 . 000000 , 2 . 000000] {(-9 .800000) / ( (p) * (-4 . 900000)+ 
(log(p))* (2.300000)+(2.300000)*(log(p))+(4.900000)*(-(p))), 
((logCp))*(2. 000000)) /((p)*(-4.900000)+(log(p))*(2. 300000)+ 

(2 . 300000) * (log(p) ) + (4 . 900000)* (- (p) ) ) , (4 . 600000) / ( (p) *(-4 . 900000) 
+ (log(p) ) * (2 . 300000)+(2 . 300000) * (log(p) ) +(4 . 900000) ♦ (- (p) ) ) , 

( (- (p) ) * (2 . 000000) ) / ( (p) ♦ (-4 . 900000) + (log (p) ) * (2 . 300000) + 

(2 . 300000) ♦ Clog(p) ) + (4 . 900000)* (- (p) ) ) } ; 

/*ADDED INPUT*/ 

p * 3.246; 


unit » mmult (){test*testinv;3-; 


/*0UTPUT*/ 

test [2 . 000000 , 2 . 000000] {3 . 246000 , 1 . 177423 , 2 . 300000,4. 900000} ; 
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test inv [2 . 000000 , 2 . 000000] {0 . 371287,-0 . 089217 , -0 . 174278 , 0 . 245959 } ; 
unit»inult [2 . 000000 , 2 . 000000] {1 . 000000 , 0 . 000000 , 0 . 000000 , 1 . 000000} ; 




Chapter 4 

Graphical Simulation 


In the previous two chapters, we saw the methodology to be followed for obtaining the 
symbolic equations of the machined surface and the implementation of the symbolic 
manipulator. But the main motive of the work is to render the surfaces graphically. 
In the following section we discuss how it is achieved. 

4.1 Surface Rendering 

There are many ways a surface can be rendered. The major ones are ; 

• Wireframe 

• Filled Polygons 

• Spline surfaces 

The first approach, though simple and fast, is not suitable for the present need. 
The surfaces are of free form, and except for very simple ones it becomes difficult to 
visualize them with this model. The last approach gives a very realistic rendering 
of the surfaces. But the calculation of splines dynamically is a very computation 
intensive operation and this fact makes them unfit for our purpose. 

The second approach, coupled with shading techniques (Goraud in the present 
case), is ideally suited for the current application. In this approach, the surface is 
approximated by polygons. In order to reduce the faceted appearance of the surface 
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and get a snioothor surface, the intensity is interpolated across polygons. For this 
purpose normals need to be specified at each vertex of each polygon in the surface. 

4.1.1 Calculation of Normals 

In general there are two approaches used for calculation of normals depending upon 
the tradeoff between the smoothness of animation and realism of the rendition. 

The first approach, fast but inaccurate, is to calculate the normals at a vertex by 
taking the cross product of two edges of the polygon incident on that vertex. The 
shading based on such normals is uneven and a close examination reveals the faceted 
nat are of the surface. 

In the second approach, the normal at a vertex is calculated by taking the resultant 
of ail normals due to all the edges incident on a vertex, instead of just two. This gives 
a much smoother shading to the surface. 

In our case, we have the choice of going for an even better rendition by calculating 
exact normals at each vertex. Since we have the symbolic equations of the cutter 
surface normals in the So frame of reference, all that is to be done is to transform 
them to the frame of reference. This is done by postmultiplying them with the 
transformation matrix Mq using symbolic manipulator. Having obtained the symbolic 
normals, we only have to give different values to u, v and t to get the normals at 
different vertices. The surfaces rendered thus are highly realistic as the photographs 
illustrate. The animation also does not suffer. 

4.2 Animation 

The field of computer animation has progressed appreciably during the past decades 
owing to evolution of special purpose hardware and powerful techniques such as back- 
face culling, hidden-surface removal and double buffering which has made it possible 
to render motion realistically. 

A manufacturing process is a dynamic process. It is not very impressive to show 
just the cutter and the machined surfaces. If the kinematics are also shown then the 
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<l«‘sjgi!(“r rail vi5ualiz<' tho motion and change the kinematics appropriately to get the 
(h'sired Kiirfare. 

Ii» this work two distinct types of animation may be observed. One is due to 
the kinematic parameters supplied earlier and is fixed. The other is because of the 
moveiiH‘t>t of the observer, which is interactive. The capability of moving the ob- 
.server enables to view the machining process from different view points. This kind of 
animation is a standard in graphics and we won’t discuss it any further than saying 
that many options are available for the user. 

'I'he first kind of animation, basically, shows the motion of the cutter as seen from 
the r<’ferrnre frame $2 and the machined surface as the blank is being cut. This is 
achieved its follows: 

I’rom the symbolic manipulator we already have the cutter surface as a function 
of u, t> and t. u and v are the parameters specifying the surface and t is the time. 
Stepping t between two instants of time with the required value of time step, we 
obtain the position of the cutter at the corresponding time as guided by the kinematic 
parameters. 

Now, to show the machined surface we compute the cutter profile satisfying the 
<"oiidition of contact and store it along with the normals. It should be clear that the 
machiiufd surface is completely calculated before any output appears on the screen. 
Only that portion of the surface is displayed which is actually cut till the present 
moment. That gives the effect of the blank being cut then and there. 

4,3 User Interface 

The user interface provided is highly customized to the present application. To use 
the symbolic manipulator for other purpose say, for computing a different algorithm, 
one has to access it separately. 

The menu is textual and the user is prompted to select one of the following four 
sub-menus (or quit): 


• Cutter specifications 
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• KiiH'inatifs 

• C»ra{>hics 

• Hun 

'I'he Cuiitr specifications menu asks for the values of nine parameters of the 
generic cutter. The Kinematics menu similarly prompts the user for the kinematics 
of the cutter as well as the blank surfaces. Graphics menu controls the rendering 
asjjects stich as the background colour , the type of normals to be used for rendering 
the cutter surface (fast or exact), the value of time step and the like. The final menu, 
Hint gives option.s for compiling, cleanup, running the program etc. 

Default values are provided in all appropriate places. 

Sir»ce it is very difficult to predict beforehand what the exact numerical values 
should be for the coefficients of various parameters, as many as 25 variable parameters 
are supported. Any variable parameter which needs to be specified at the run time 
should have a name of the type varparamX , where X is any number between 1 and 
25, both included. This limit of 25 on variable parameters is purely arbitrary and can 
be set to any desirable value by minimal addition of code. These parameters should 
be used sparingly as it would increase demands on system memory and time. 


4.4 Examples 

In this section we give three examples, which have been directly tested on HP-TSRX 
workstations. The output (graphical) is shown in the photographs that follow: 
Example:!: 


Is arc 

tangential? 

d 

r 

a 

b 

hi 

e 

f 

h 

Yes 

5.181 

4.0 

0.0 

24.141 

20.0 

0.0 

4.0 

50.0 


Table 4.1: The Cutter Parameters for the Conical Cutter 


Afo. 
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}-xi>lanati{>n: The cutter is stationary with its axis coinciding with the Y-axis of 
the global C(K)rdinate syvstem. The blank is rotating around X-axis of global coordinate 
system and at the same time moving along the positive X-axis (of ^o). See Tables(4.1, 
4.2 and 4.3). 

Exampk:2: 

Xi(tX, u) = ucos v 
yi(u,v) = usinu 
zi(u,v) = au^ -i- bu + c 


Where 

0 <u < 5.339036 + 01 

0 < V < 27r 

a = 1.68979e-02 
6 = 1.02956e-06 
c = 4.91379e + 01 

I'he Cutter Equations for End Milling Cutter 

Explanation: The cutter axis is in the XY plane, making an angle of ^ with the 
positive Y-axis (of Sq frame of reference) and the tip of the cutter is at a distance of 
2.0 units along Y-axis initially. 

The cutter is moving away in positive Y direction while maintaining the above 
configuration. See Table(4.4). 

The blank is rotating around X-axis of global coordinate system and at the same 
time moving along the positive X-axis (of So)- See Table(4.5). 

Example:3: 

Explanation: Initially the cutter tip is at a distance of 200.0 units from the origin 
along positive Y direction and moves sinusoidally along the Y-axis. See Table(4.6). 

The blank rotates around X-axis of global coordinate system and at the same time 
moves along the positive X-axis (of So)- See Tables(4.7 and 4.8). 
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h 

Cl 


<f>i 


0.0 

0.0 

0.0 

■■- 2.0 

0.0 

0.0 


I able 4.2: The Cutter Kinematics of the Conical Cutter 


02 


C2 

$2 

<f >2 

^2 

-50.O7-.i»Q,O125 

— 2>Qfg: 

0.0 

0.0 

t*0.0l25 

0.0 

0.0 


Table 4.3: The Blank Kinematics for the Conical Cutter 


Cl 

bi 

Cl 


4^1 

01 

0.0 

2.0+0. ll*t 

0.0 

TT 

7.0 

0.0 

TT 

12.Q 


I able 4.4: The Cutter Kinematics of the End Milling Cutter 


02 

62 

^2 

^2 

02 

02 






0.0 


Table 4.5: The Blank Kinematics for the End Milling Cutter 


Is arc 

tangential? 

d 

r 

a 

b 

h. 

e 

f 

h 

'iV.s 

25T 

4.0 

-10.0 

-5.0 

50.0 

8.061 

2.64 

20.0 


Table 4.6: The Cutter Parameters for Bottom Angle Cutter 


Cl 

by 

Cl 

1 1 


01 

0.0 

100.0*(1.0+cos(4.0*t*j^)) 

0.0 


0.0 

0.0 


Table 4.7: The Cutter Kinematics of the Bottom Angle Cutter 


02 

62 

C2 

62 

02 

02 




t*0.0125 

0.0 

0.0 


Table 4.8: The Blank Kinematics for the Bottom Angle Cutter 



























Figure 4.2: The Helical Screw being Machined 
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Figure 4.5; The Spiral Screw being Machined 



Figure 4.6: The Spiral Screw Surface 
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Figure 4.7: The Machining of Star Shaped Surface 



Figure 4.8: The Star Shaped Surface 





Chapter 5 
Conclusions 


The present work encompasses a fairly large number of cases encountered in practice. 
It renders correctly all such surfaces which axe cut by axisymmetric cutters such 
that a portion once cut is not cut again. Even cases in which the cutters are not 
axisymmetric, can be handled by taking their biparametric equations directly. The 
cases where a satisfactory rendering could not be obtained are considered in the 
following section. 

5.1 Considerations for Further Work 

I'he main problem with the present approach crops up because of the self intersection 
of the cutter path. Such cases are not very rare in practice. These cases have to be 
solved computationally, as it may not always be feasible to solve them symbolically. 

To do this one considers the cutter surfaces at a very close interval of time. The 
snapshots so obtained are checked for intersection and the portions which are inside 
are removed. The resultant surface is still not in an acceptable form as even for 
very close snapshots, the machined surface obtained by the above method might 
show severe irregularities not to be seen in practice. Such irregularities can then be 
smoothed out using interpolation near the points of intersection of surfaces. 

It may be noted that the intersection need not be with two successive snapshots 
only. It can be with the machined surface cut long before. Therefore in a general 
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ra5f‘, one of the surfaces (the machined surface, to be precise), will be an approximate 
surface based on the points stored for rendition. 

All through the work we have assumed that the cutter is fully immersed in the 
blank which is generally not the case. A further extension could therefore be to define 
an initial shape for the blank. Then to render the machined surface, we proceed as 
before and at the end check the developed surface for intersection with the blank 
surface. The portions which are outside are then removed leaving the actual machined 
surface. 


5.2 Some Techniques and Improvements 

In this section we will discuss how to add more functionality to the symbolic manip- 
ulator and how to make it more efficient in terms of memory and minimization of 
expressions. 

As promised earlier, we will first look into how substitution of symbols can be 
made partial. The function substitute first substitutes a symbol and then substitutes 
the previously substituted expression. All that is to be done is to replace this code 
with an iterative loop whose limit will specify how deep the substitution should be. 

I’he next issue is how to add more functionality to the symbolic manipulator. 
The manipulation of functions is mainly localized in three routines, namely reserved, 
which checks whether a function is a procedure supplied with the symbolic manipu- 
lator, evaluate, which on every occurence of a function name, calls the corresponding 
module in the program and last of all Display (and its clones Sisplay and /display), 
which formats the output. To add a function one writes the function code and: 

• Puts the name in function reserved 

• Translates the function usage to function call in evaluate 

• Formats the output as needed in Display (and other clones) 
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Improvements: The symbolic manipulator is very efficient in terms of memory 
and minimization. But still, in future this might not be adequate. There is scope for 
further improvement in these areas. 

1) Memory: In the design it was specified that (except for inverse) the symbolic 
manipulator frees the Table only between two declarations. It is possible to free some 
portion of the Table even inside a single declaration. For this purpose we can either 
reuse or free the Tokens which are no longer a part of the expression. This has to be 
done carefully because even if a Token is currently not a part of the expression tree, 
it may be later. In such cases one can lock those Tokens by making the valid field, 
say, 1 and releasing the other unused Tokens. 

2) Minimization: The present implementation recognizes only 2” of n! combina- 
tions as identical in case of commutative operators. To recognize all n! combinations, 
one can sort the expressions on the key provided in the key field of the listnode and 
then compare the expressions. For the cases of factors, it is possible to de-factorize 
the expressions, sort them and compare. These two improvements are sure to increase 
its memory efficiency too, as the two are interrelated. In addition there is hardly any 
bound on the trigonometric and other simplifications possible. 
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Appendix A 


Formal Specification of Symbolic 
Manipulator 


The following is the formal specification of the symbolic manipulator. 


%% 


Xtoken 

IDENTIFIER 

Xtoken 

CONSTANT 

Xtoken 

EQ_0P NE.OP 

Xtoken 

SIN COS TAN COT SEC COSEC ASIN ACOS ATAN 

% token 

ACOT ASEC ACOSEC SINK COSH 

Xtoken 

TANK COTH SECH COSECH ASINH ACOSH ATANH 

%token 

ACOTH ASECH ACOSECH 

% token 

SORT MATRIX LOG ABS 

Xstart 

bgn 

n 



primary.expr 


: identifier 
I CONSTANT 
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expo.expr 


©xpo_operator 


unary.expr 


unary .operator 


' ( ' expr ' ) ' 

function.def inition 


primary.expr 

primary.expr expo.operator expo.expr 


> 


> 


unary .operator cast .expr 


SIN 

COS 

TAN 

COT 

SEC 

COSEC 

ASIN 

ACOS 

ATAN 

ACOT 

ASEC 

ACOSEC 

SINK 

COSH 

TANK 


COTH 


SECH 


cast_expr 


multiplicative 


additive.expr 


COSECH 
AS INK 
ACOSH 
ATANH 
ACOTH 
ASECH 
ACOSECH 
LOG 
SQRT 
ABS 
> + > 


expo_eipr 

unaxy_expr 


,expr 


cast_expr 

multiplicative_expr 
multiplicative.expr ’/' 
multiplicative.expr ’%' 


cast_expr 

cast.expr 

cast.expr 


: multiplicative.expr 

1 additive.expr multiplicative.expr 

I additive.expr multiplicative.expr 
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equal ity_expr : 

I 

1 


assignment_expr : 

I 


assigiment_operator : 


expr 

declaration 


init_declarator 


inat_init_declarator 


additive.expr 

equal ity_expr EQ_DP additive_expr 
equal ity.expr NE_DP additive.expr 


equal ity_expr 

unary_expr assignment_operator assignment.expr 


J 


as s ignment _ expr 


: iny_type_spec mat.init ..declarator ’ ; ’ 
I init .declarator ’ ; ’ 


: declarator 

I declarator initializer 

I 

: mat.declarator 

I mat.declarator mat. initializer 


declarator 


aat_declaxator 


By_const_list 


identif ier_list 


mat. initializer 


mat_initializer. 


initializer 


: identifier 

I declarator '(’ identif ier_list ')' 
I declarator '(' 


: identifier ’ [’ my_const_list ']’ 


: CONSTANT 

I my_const_list ' , ’ CONSTANT 


: identifier 

I identif ier_li St ' , ’ identifier 


: expr 

I mat_initializer_list 


list: mat_initializer 

I mat.initializer.list mat.initializer 


: expr 

I initializer '}' 


compound.statement : expression.statement ’}' 
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expression.statement : ' ; ' 

I expr ’ ; ' 


bgn ; file 

I 

: external_def inition 
I file external_def inition 


eiternal_def inition : function.definition 

I declaration 


function.definition : declarator function_body 


my_type_spec : MATRIX 


function_body 


: compound.statement 


identifier : IDENTIFIER 


n 



Appendix B 


A Sample Input to Symbolic 
Manipulator 


Here wc give a sample input of the symbolic manipulator. 

PI - 3.14159265; 

DEG- PI/180.0; 

PZERO - l.OE-10; 

INF « 1.0E+20: 

TRUE - 1.0; 

FALSE - 0.0; 

/♦CUTTER PARAMETERS*/ 

IFR - TRUE; 
d - 1.0; 
r « PZERO; 
a»PZER0; 
b-0.0; 
lil“70.0; 
ein.it -0.5; 
finit * 0.0; 
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hinit«0.0; 

/♦KINEMATICS*/ 

thetal - -PI/2.0; 

phil ■ 0.0; 

xil » PI/varparam8; 

tlieta2 ■ t * varparaml ; 
plii2 » 0.0; 
xi2 » 0.0; 

al»0.0; 

bl* varparamS ♦ (1.0 + varparam5*t + cos(varpaxafli4*t*PI/180.0)) ; 
cl»0.0; 

a2“ -varparain2 ♦ t * varparaml/ (2.0 ♦ PI); 

b2-0.0; 

c2«0.0; 


boolb » (abs(b)/(abs(b) + 1 .OE-20) ) *0.00001 ; 
noth « abs (boolb - 1.0); 

e*einit ♦ abs(IFR - 1.0) + IFR* (d/2.0 - r*((cos(a) - 
8in(b))/cos(a + b))); 

f«finit * abs (IFR - 1.0) + IFR*(r + e*siii(a))/cos(a) ; 
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iPc- 

(o + f*tan(a))*cos(a)*cos(a) - sqrt(abs(((e + f ♦tamCa) )*cos(a) 
*cos(a)) "2 .0 - (e*e + f*f - r*r)*cos(a)*cos(a)) ) ; 

zPc*xPc*t2ui(a) ; 


zQc » 

(f + (e + (d/2.0)*(tan(a)*tan(b) - 1.0))*t2m(b))*cos(b)*cos(b) 

+ sqrt(abs(((f + (e + (d/2.0)*(taii(a)*taii(b) - 1 .0))*taii(b))*cos(b) 
♦cos(b))“2.0 - (f*f + (e + (d/2.0)*(taji(a)*taii(b) - 1.0))'‘2.0 - r*r) 
*cos(b)*cos(b))) ; 

xQc»zQc*tan(b) - (d/2.0)*(tan(a)*taii(b) - 1.0); 

b « IFR+aotb+zQc + abs(IFR*notb - 1.0)*hiiiit; 


xSc»linaii(b) - (d/2.0)*taii(a)*tan(b) + (d/2.0); 

zSc*h. ; 


8l»=sqrt(abs(iPc*xPc + zPc*zPc)); 

matrix tantempl [1 .0,2.0] ={zQc - f , xQc -e}; 


matrix tantemp2[1.0,2.0]={zPc - f , xPc -e}; 
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s2“ab8(r*( taninv(){tantempl ;} - taninv(){tzmtemp2;})) ; 

83*8qrt(ab8((xSc - xQc)“2.0 + (zSc - zQc)‘2.0)); 
s4»'abs(hl) ; 

matrix cPul [1.0,4.0]*-Cxcl,ycl,zcl,l .0>; 
matrix cPu2 [1.0,4.03 ={xc2,yc2 , zc2 , 1 . 0> ; 
matrix cPu3[1.0,4.0]={ic3,yc3,zc3,l .0}; 
matrix cPii4[1.0,4.0]*-Cxc4,yc4,zc4,l .0}; 

xcl*u*cos(a) ; 

ycl»0.0; 

zcl*u*siii(a) ; 

matrix taiit6mp3[1.0,2.0]*{sl*sin(a) - f , sl*cos(a) - e}; 
xc2«e + abs(r)*cos(«.taiiiiiv(){tajiteinp3;} + ((u - sl)/r)); 
yc2»0 . 0 ; 


zc2*f + abs(r)*sin(__taiiinv(){tantemp3;} + ((u - sl)/r)); 




xc3“(u - (si + s2 + s3))*sin(b) + iSc ; 

yc3*0.0; 

zc3»(u - (si + 82 + s3))*cos(b) + zSc ; 

xc4*xSc; 

yc4“0 . 0 ; 

zc4*u “ (si + s2 + s3) + zSc ; 
matrix MlcC4.0,4.0]* 

{cos (v) .sin (v) ,0 .0 ,0 .0,-sin (v) ,cos(v) ,0.0, 
0.0, 0.0, 0.0, 1.0, 0.0, 0.0. 0.0, 0.0, 1.0}; 

pll» imnult(){cPul*Mlc ;}; 

pi 2* miault(){cPu2*Mlc ;}; 

pl3» mmult (){ cPti3*M1c ;}; 

pl4* iiunult(){cPu4*Mlc ;}; 

matrix pll [1 • 0 ,4 . 03 *{u*cos(v) ,u*sin(v) , a*u*u 


b^u + 0,1.0} 


matrix pl2[l .0,4.0]*{u*cos(v) ,u*sin(v) .a*n*u + b*u + c.l.O} 
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matrix pl3[l . 0 ,4.0]»{u*cos(v) ,u*sin(v) , a*u*u + b*u + c,1.0}-; 
matrix pl4[l . 0,4.0] »{u*co8(v) ,u*siii(v) ,a*\i*u + b*u + c,1.0]-; 
itheta2 * -tlieta2 ; 


iphi2 » -phi2; 

ixi2 ■ -xi2; 


ia2« -a2; 

ib2« -b2; 

ic2» -c2; 

matrix R01x[4.0,4.0]« 

{1.0,0.0,0.0,0.0,0.0,cos(thetal) .sinCthetal) 

,0.0, 0.0, -sin. (tbetal ) , cos (thetal ), 0.0, 0.0, 0.0, 0.0, 1.0}; 

matrix R01y[4.0,4.0]= 

{cos(pbil) ,0.0,-sin(pbil) ,0.0,0. 0,1. 0,0.0, 

O.O.sinCpbil) ,0.0,cos(phil) ,0.0,0.0,0.0,0.0,1.0}; 
matrix R01zC4.0,4.0]* 

{co8(xil),sin(xil),0.0,0.0,-sin(xil),cos(xil) 

, 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 }; 


matrix TOl [4.0,4. 0] 




— — ^ 

{ 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 

0.0.al,bl,cl,1.0>; 

MOl • — minult(){R01x*__mmult(){R01y*._minult(){R01z*T01;};};}; 

HOld * differentiate(t){M01;}; 

matrix R02x[4.0,4.0]= 

{1 . 0, 0. 0, 0. 0,0. 0,0.0, cos (theta2) ,sin(theta2) , 
0.0,0.0,-sin(theta2),cos(theta2) ,0.0, 0.0, 0.0, 0.0, 1.0}; 

matrix R02y [4.0,4. 0]= 

{cos (plii2) , 0 . 0 , -sin (p]ii2) , 0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 0 . 0 , 
sin(phi2) ,0 .0,cos(phi2) ,0.0, 0.0, 0.0, 0.0, 1.0}; 

matrix R02z[4.0,4.0]* 

{cos (xi2) , 8in(xi2) ,0.0,0. 0 ,-sin(xi2) , cos(ii2) , 

0.0, 0.0, 0,0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 1.0}; 

matrix T02[4.0,4.0]» 

{ 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 

0.0,a2,b2,c2,1.0}; 

M02 » __nimult(){R02x* mmult(){R02y*__nmult(){R02z*T02;};};}; 

M02d » differentiate(t){M02;}; 

matrix iR02x[4. 0,4.0]“ 

{1.0, 0.0, 0.0, 0.0, 0.0, cos(itlieta2) , sin(itlieta2) 

,0.0,0.0,-sin(itheta2) ,cos(itheta2) ,0.0, 0.0, 0.0, 0.0, 1.0}; 
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matrix iR02yC4.0,4.0]» 

{co8(iphi2) ,0.0,-sin(iphi2),0.0,0.0.1.0,0.0, 

0.0,sin(iphi2) ,0 .O.cos(iphi2) ,0.0,0 .0,0. 0,0.0, 1 .0}; 

matrix iR02z[4.0,4.0]= 

{cos(ixi2) ,siii(ixi2) ,0.0,0.0,-sin(ixi2) ,cos(iii2) 

, 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 >; 

matrix iT02[4.0,4.0]» 

{ 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 

ia2,ib2,ic2,1.0}; 

M20 « __iiiverse(){M02;}; 

/♦M20 » mnnilt(){iT02* minult () {iR02z* ininult(){iR02y*iR02x 

would be more efficieut in terms of memory*/ 

p21tmp* _jmmult(){M01*M20 ;>; 

p21 * mmult(){pll*p21tmp;>; 

p22 * mmult(){pl2*p21tmp;>; 

p23 * mmult (){pl3*p21tmp;} ; 

p24 « __minult(){pl4*p21tmp;>; 


Nlpll » cross 


(){__dif f erentiate (u) {pll ;}*__dif f erent iate (v) {pi 1 ; } ; > ; 
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Nlpl2 * _ .cross (){ — differentiate(u){pl2;}* differentiate(v)-Cpl2;};3'; 

Nlpl3 » — cross (){ differeiitiate(u){pl3;}* differentiate(v)-Cpl3;3-;3-; 

Nlpl4 * cro8s()-C differeiitiate(u){pl4;}* differentiate(v)-Cpl4;} ;}; 

NOpll»__in2Dult(){Nlpll*M01;>; 

N0pl2»_.nmiult(){Nlpl2*M01;}; 

N0pl3»._minult ( ) {Nlpl3*M0 1 ; } ; 

N0pl4»„mmult (){Nlpl4*M01 ; }; 

N2pl 1* mmult ( ) {NOpl 1*M20 ; } ; 

H2pl2-..Bmult 0 {N0pl2*M20 ; } ; 

N2pl3»__ininult (){N0pl3*M20 ; >; 

N2pl4*' mmult () {N0pl4*M20 ; }■ ; 

Vlltmpm ■ mmult ()-Cpll*M01d;}; 

V21tmpm » mmult()-Cp21*M02d;}; 

V0pl21 »__msub()-CV21tmpm - Vlltmpm;}; 


V12tmpm * mmult ()'Cpl2*M01d; } ; 
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V22tiBpB - __Bunult(){p22*M02d;}; 

V0pl22 ■ — msub()’{V22tffl[pm - V12tmpm;>; 

ViStBpB ■ Bmult()-Cpl3»M01(i;}; 

V23tapia « iHinult(){p23*M02d;}; 

V0pl23 “_^ttsub()<V23tmpm - V13tmpm;>; 

V14t«pm * nmTilt(){pl4*M01d;}; 

V24tmpm * __mmTilt()-Cp24*M02d;}; 

V0pl24 » msub(){V24tmpm - V14tmpm;>; 

cod* dot(){N0pll>*‘V0pl21;}; 

coc2» ._dot(){N0pl2*V0pl22;}; 
coc3* dot(){N0pl3*V0pl23;}; 


coc4» dot(){N0pl4*V0pl24;}; 



